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FOREWORD 

This  work  represents  part  of  six  man-months  of  effort,  the  ultimate  objective 
of  which  is  the  direct  prediction  of  the  onset  and  details  of  submarine  pressure 
hull  rupture  caused  by  underwater  explosive  attack.  The  immediate  objective  of 
this  work  is  to  consider  an  approximation  offering  some  hope  of  circumventing 
the  difficulties  of  applying  the  Doubly  Asymptotic  Approximation  (DAA),  an 
important  but  limited  analytical  tool  for  treating  fluid-structure  interaction  in 
cases  involving  acoustic  shock  waves  impinging  on  underwater  structures. 
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INTRODUCTION 

The  Doubly  Asymptotic  Approximation,  as  Implemented  In  USA-STAGS,  has  been 
successfully  used  to  model  fluid-structure  interaction  in  a  number  of  studies 
of  the  response  of  underwater  structures  to  explosions.  Unfortunately  the  USA 
implementation  of  the  DAA  appears  to  embody  an  utterly  unrealistic  feature 
according  to  which  negative  fluid  pressures  arise  and  exert  a  kind  of  suction 
force  on  the  structure.  In  fact,  the  DAA  should  only  be  applied  up  to  the  time 
at  which  the  total  pressure  becomes  zero.  Furthermore,  this  time  occurs  extremely 
rapidly  in  typical  explosion  loading  situations. 

In  this  report  we  consider  a  simple  approximation  to  the  DAA,  which  takes 
advantage  of  the  rapid  pressure  cutoff. 

MODEL  PROBLEM 

We  consider  a  plane  acoustic  wave  impinging  on  an  infinite  unrestrained 
air-backed  plate  of  mass  m  per  unit  area,  as  illustrated  in  Figure  1.  The 
incident  pressure  is  subject  to  decay  typical  of  underwater  explosions: 

Pi  =  P0  exp  (-t/e). 

The  total  pressure  on  the  plate  is  the  sum  of  the  incident  pressure,  the  reflected 
pressure,  and  the  relief  pressure  arising  from  plate  motion: 

PT  =  Pi  +  Pr  +  Pr-  (1) 

The  DAA  uses  tlv.  approximation  [1] 

PR  =  Pr  Pr  =  -pcv 
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where  r  is  the  water  density,  c  Is  the  acoustic  wave  speed  In  water,  and  v  Is 
the  plate  velocity.  The  equilibrium  equation  becomes 

dv 

mdt  x  Pt 

*  2p0  exp(-t/o)  -pcv  (2) 

Now  using  the  Jlmensionless  quantities 

w  *  v/c  t  *  t/e 

k  =  pce/m  v  =  2pQ9/cm 

equation  (2)  is  rewritten  as 

+  <w  3  v  exp(-t)  (3) 

Equation  (1)  is  called  the  Doubly  Asymptotic  Approximation  since  it  is 
expected  to  give  correct  answers  for  very  short  and  very  long  times. 

The  further  simplification  to  be  studied  in  this  report  is  to  neglect 
relief  pressure  except  for  determining  the  cutoff  time  and  velocity.  The 
rationale  underlying  this  approximation  is  illustrated  in  Figure  2.  Owing  to 
pi  ate  inertia,  appreciable  plate  motion  and  therefore  appreciable  relief  do  not 
occur  until  some  time  has  passed.  But  by  this  time  the  major  part  of  the  pressure 
decay  is  likely  to  have  occurred.  If  so,  the  cutoff  pressures  <wr  calculated 
using  relief  and  kw^  without  relief  both  intersect  the  curve  for  incident  plus 
reflected  pressure  in  tho  flat  part  of  this  curve.  In  this  event,  the  difference 
between  the  intersection  values  of  icwr  and  <w^  should  be  small. 

CASE  1:  RELIEF  PRESSURE  RETAINED 

With  the  initial  condition  w(o)  =  o,  the  solution  of  (3)  is 

w  =  — —  [exp(-icT)  -  exp(-x)]  (4) 
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Cutoff  occurs  when 


pcv  *  2p0  exp(-t/c) 


and  hence  when 

kw  *  exp(-x).  (5) 

Substituting  (5)  and  (4)  and  solving  for  the  cutoff  time  tc  leads  to 

tc  a  -  tnU1'*)  (6) 

The  cutoff  velocity  wr  is  now 

wr  *  ~  exp(-r) 

(dir) 

*  \>  K 

CASE  2:  NEGLECT  RELIEF  PRESSURE  EXCEPT  AT  CUTOFF 

In  this  approximate  treatment  the  relief  pressure  in  the  governing  equation 
is  neglected,  leading  to 

»  v  exp(-t) 

with  the  solution 

w  =  v(l-exp(-t) 

Now  cutoff  still  occurs  when  the  actual  relief  pressure  causes  the  total  pressure 
to  vanish; 

pcv  =  2p0  exp(-t/e) 

giving 

<w  =  v  exp(-t) 

The  cutoff  time  and  velocity  are  now 

xc  =  -  8,n  (</ ( 1  +  <)) 
wN  =  v/(l  +  e) 

The  cutoff  times  and  velocities  are  illustrated  in  Figure  2  for  the  cases 
with  and  without  relief. 
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COMPARISON  OF  ACTUAL  AND  APPROXIMATE  CUTOFFS 

To  measure  the  error  of  neglecting  relief,  we  form  the  ratio  wr/w^: 

(JL-) 

Wr/WN  «  (1  +  If)*1'*  3  <fr(ic) 

The  limiting  behavior  of  $(*:)  is  as  follows: 

lim  4>(*c )  »  1 

K  0 

lim  $(k)  *  1 

lim  a(K)  =  2/e  *  .736 

K  -*■  1 

Also,  the  minimum  value  of  $(<)  is  <j>(l)  =..736. 

Figure  3  gives  computational  results  for  the  dependence  of  $(<)  on  k. 

• 

EXAMPLES 

As  in  reference  2  we  consider  two  one  inch  thick  target  plates,  one  made 
of  steel  and  the  second  made  of  titanium.  We  also  consider  two  warhead  types; 
the  first  involves  50  lbs  of  TNT  at  300  feet  standoff  [3]  while  the  second 
consists  of  a  copper  lined  shaped  charge  fired  under  water,  on  which  we  found 
some  data  [2].  The  values  of  k  and  $(<)  can  be  organized  into  the  following  table. 


Conventional 

Copper  Jet 

K.  ${<) 

<»  $(<) 

Steel 

5.38,  0.31 

0.27,  0.78 

Titanium 

8.97,  0.84 

0,45,  0.75 

As  Figure  3  makes  clear,  the  maximum  error  in  cutoff  velocity  cannot  exceed 
25%.  Furthermore,  the  error  is  approximately  twenty  percent  in  the  range  of 
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Interest.  It  therefore  appears  possible  to  assume  that  the  true  cutoff  velocity 
vc  is  approximately  .8v^  where  v^  is  the  cutoff  velocity  calculated  assuming 
no  relief. 

DISCUSSION 

We  have  used  a  model  problem  to  examine  the  error  of  a  simple  approximation 
obviating  the  calculation  of  the  relief  pressures  used  in  the  Doubly  Asymptotic 
Approximation.  The  error  is  a  function  of  m/pce,  which  may  be  interpreted  as 
the  ratio  of  structural  response  times  to  fluid  pressure  decay  times.  The 
greatest  relative  error  is  about  27%  and  it  occurs  when  m/pce  =  1.  For  practical 
examples  involving  steel,  titanium,  conventional  warheads,  and  shaped  charge 
warheads,  it  appears  that  the  error  is  typically  about  20%,  suggesting  a 
"correction  factor." 

The  modest  error  of  the  present  approximation  in  respect  to  cutoff  velocity 
is  of  interest  in  its  own  right,  particularly  in  the  common  case  in  which  most 
of  the  energy  available  from  the  shock  wave  is  transferred  to  the  tarqet  plate 
before  the  true  cutoff  time.  However,  some  energy  is  transferred  to  the  target 
plate  after  cutoff  by  the  fluid  still  flowing  toward  the  plate.  We  hope  to 
obtain  a  simple  approximation  for  this  part  of  the  process  in  our  future 
investigations. 
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FIGURE  1  TARGET  CONFIGURATION 
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CUTOFF  VELOCITY  WITH  RELIEF  TO  CUTOFF  VELOCITY  WITHOUT  RELIEF 
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